Pygmy dipole resonance: collective features and symmetry energy effects 
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A very important open question related to the pygmy dipole resonance is about its quite elusive 
collective nature. In this paper, within a harmonic oscillator shell model, generalizing an approach 
introduced by Brink, we first identify the dipole normal modes in neutron rich nuclei and derive the 
energy weighted sum rule exhausted by the pygmy dipole resonance. Then solving numerically the 
self-consistent Landau- Vlasov kinetic equations for neutrons and protons with specific initial condi- 
tions, we explore the structure of the different dipole vibrations in the Sn system and investigate 
their dependence on the symmetry energy. We evidence the existence of a distinctive collective 
isoscalar-like mode with an energy well below the Giant Dipole Resonance (GDR), very weakly 
dependent on the isovector part of the nuclear effective interaction. At variance the corresponding 
strength is rather sensitive to the behavior of the symmetry energy below saturation, which rules 
the number of excess neutrons in the nuclear surface. 
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One of the important tasks in many-body physics is 
to understand the emergence of the collective features as 
well as their structure in terms of the individual motion 
of the constituents. The steady progress of experimen- 
tal methods of investigation opens now the possibility 
to study very neutron rich nuclei, beyond the limits of 
stability. The goal is to have a unified picture of the evo- 
lution of various nuclear properties with mass and isospin 
and to test the validity of our theoretical understanding 
over an extended domain of analysis. 

New exotic collective excitations show up when one 
moves away from the valley of stability [1]. Their ex- 
perimental characterization and theoretical description 
is a challenge for modern nuclear physics. Recent exper- 
iments provided several evidences about their existence 
but the available information is still incomplete and their 
nature is a matter of debate. 

An interesting exotic mode is the Pygmy Dipole Res- 
onance (PDR) which was observed as an unusually large 
concentration of the dipole response at energies clearly 
below the values associated with the GDR. The latter 
is one of the most prominent and robust collective mo- 
tions, present in all nuclei, whose centroid position varies, 
for medium-heavy nuclei, as SQA~^/^MeV . Adrich et 
al. [2!] reported the observation of a resonant-like shape 
distribution with a pronounced peak around IQA'IeV in 
^^^Sn and ^^"^Sn isotopes. A concentration of dipole ex- 
citations near and below the particle emission thresh- 
old was also observed in stable Sn nuclei, a systemat- 
ics of PDR in these systems being presented in It 
was concluded that the strongest transitions locate at 
energies between 5 and ^.bMeV and a sizable fraction of 
the Energy- Weighted Sum Rule (EWSR) is exhausted by 
these states. From a comparison of the available data for 
stable and unstable Sn isotopes a correlation between the 
fraction of pygmy strength and isospin asymmetry was 



noted In general the exhausted sum-rule increases 
with the proton-to-neutron asymmetry. This behavior 
was related to the symmetry energy properties below sat- 
uration and therefore connected to the size of the neutron 
skin 0-01 • However other theoretical analyses suggest a 
weak connection between the PDR and skin thickness ^] . 

In spite of the theoretical progress in the interpreta- 
tion of this mode within phenomenological studies based 
on hydrodynamical equations 0, , non-relativistic mi- 
croscopic approaches using Random Phase Approxima- 
tion (RPA) with various effective interactions [ll| - [l3| or 
relativistic quasi-particle RPA [13, [13 j and new exper- 
imental information p^| - [T9| . a number of critical ques- 
tions concerning the nature of the PDR still remains. 
This includes the macroscopic picture of neutron and pro- 
ton vibrations, the exact location of the PDR excitation 
energy and the degree of collectivity of the low-energy 
dipole states, the role of the symmetry energy [lO]- Some 
microscopic studies predict a large fragmentation of the 
GDR strength and the absence of collective states in the 
low-lying region in ^^'^Sn [2l| . 

The purpose of this letter is to address the impor- 
tant issue related to the collective nature of PDR. In the 
first part, within the Harmonic Oscillator Shell Model 
(HOSM) for neutron rich nuclei, we show that the co- 
ordinates associated with the neutron excess vibration 
against the core, and to the dipole core mode respec- 
tively, are separable and derive the EWSR exhausted by 
each of them. Then we adopt a description based on the 
Fermi liquid theory with effective interactions and inves- 
tigate the dynamics and the interplay between the dipole 
modes identified in HOSM. This self-consistent treatment 
allows us to inquire on the role of the symmetry energy 
and its density dependence upon the dipole response. 

In a seminal paper , Brink has shown that for a sys- 
tem oi A = N + Z nucleons moving in a harmonic oscilla- 
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it is possible to perform a separation in four indepen- 
dent parts = Hn int + Hp mt + HcM + Ho ■ The 
first two terms determine the internal motion of protons 
and neutrons respectively, depending only on proton- 
proton and neutron-neutron relative coordinates. The 
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Hamiltonian Hcm = —. — Pcm + —^^Rni^r charac- 
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terizes the nucleus center of mass (CM) motion, while 
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X , describes a Goldhaber- 



2raNZ 2A 
Teller (G-T) ^] vibration of protons against neutrons. 

The oscillator constant K can be determined by fitting 
the nuclear size In the expressions above X and 

RcM denote the neutron-proton relative coordinate and 
the center of mass position, and we have introduced 
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the conjugate momenta P = — i-y—Pz Pn) and 
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Pcm = Pz+Pn, where Pz {Pn) are proton (neutron) to- 
tal momenta. Correspondingly, the eigenstates of the nu- 
cleus are represented as a product of four wave functions 
* = V'n intXp inta{RcM)P{X), simultaucous eigenvec- 
tors of the four Hamiltonians constructed above. For an 
El absorption a G-T collective motion with a specific lin- 
ear combination of single particle excitations is produced 
and the wave function P{X) is changing from the ground 
state to the one GDR phonon state. Denoting by Ei the 
energy eigenvalues of the system and by D the dipole 
operator , with the help of the Thomas-Runke-Kuhn 
(TRK) sum rule, the total absorption cross section is 
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Now let us turn to the physical situation, correspond- 
ing to the case of very neutron-rich nuclei, where the 
system is conveniently described in terms of a bound 
core containing all protons and Nc neutrons, plus some 
(less bound) excess neutrons N^. Thus the total neu- 
tron number N is split into the sum N = Nc + Ne 
and we denote hy A^ = Z + the number of nu- 
cleons contained in the core. In this case we have 
worked out an exact separation of the HOSM Hamilto- 
nian in a sum of six independent (commuting) quantities: 

Hsm = int + Hp^ int + int + HcM + He + Hy. 

The first three terms contain only relative coordinates 
and momenta among nucleons of each ensemble, i.e. core 
neutrons, core protons and excess neutrons and, as be- 
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fore, describe their internal motion. He = .'i — P? ' 
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X^, characterizes the core dipole vibration, while 



the relative motion of the excess neutrons against the 
core, usually associated with the pygmy mode, is de- 



termined by Hy 
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denotes the distance between neutron and proton cen- 
ters of mass in the core, while Y is the distance be- 
tween the core center of mass and the center of mass 
of the excess neutrons. The corresponding canonically 
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conjugate momenta are Pc — — — [ — Pz — ^ "^c)i and 
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of the EWSR is ex- 



Pjv^). The eigenstates 

of He and Hy are describing two independent collec- 
tive excitations and both of them will contribute to the 
dipole response since the total dipole momentum can 

be expressed as: D ^ ^X ^ ^X. + . 

Dc + Dy In this picture the PDR results in a collec- 
tive motion of G-T type with the excess neutrons oscil- 
lating against the core. The E\ absorption leads also 
to the change of the wave function associated with the 
coordinate F. ^ The total cross section for the PDR 

is: cjy = ^^i(0|[i?,y,[i/,™,Aj]|0) = ^cTD. This 
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shows that a fraction L — — 

^ NAc 

hausted by the pygmy mode. It is worth to mention that 
this result is consistent with the molecular sum rule intro- 
duced by Alhassid et al. For the tin isotope ^^"^Sn, 
if the excess neutrons were simply defined as the differ- 
ence between neutron and proton numbers, i.e. Nf. = 32, 
one would expect fy — 19.5%. This is greater than the 
value estimated experimentally, which is around 5%. A 
possible explanation for this difference is that only a part 
of the excess neutrons, A^j^, with Ny < Ne, contribute to 
PDR, the rest being still bound to the core [26| . 

Therefore it is important to test this assumption within 
a more sophisticated analysis of the dipole response. In- 
deed, a more accurate picture of the GDR in nuclei cor- 
responds to an admixture of G-T and Stenweidel- Jensen 
(S-J) vibrations. The latter, in symmetric nuclear mat- 
ter, is a volume type oscillation of the isovector density 
Pi = Pn—Pp keeping the total density p = Pn+Pp constant 
[l^j. A microscopic, self-consistent study of the collective 
features and of the role of the nuclear effective interac- 
tion upon the PDR can be performed within the Landau 
theory of Fermi liquids. This is based on two coupled 
Landau- Vlasov kinetic equations for neutron and proton 
one-body distribution functions ,fq{f,p, t) with q ^ n,p: 
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dt m dr 
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dr (?p 



(1) 



and was applied quite successfully in describing various 
features of the GDR, including pre-equilibrium dipole ex- 



citation in fusion reactions [28| . Within a linear response 
approach, it was also considered to investigate proper- 
ties of the PDR 29]. However, it should be noticed that 
within such a semi-classical description shell effects are 
absent, certainly important in shaping the fine structure 
of the dipole response |33|. We neglect here the two- 
body collisions effects and hence the main ingredient of 
the dynamics is the nuclear mean-field, for which we 



3 




(fm) 



FIG. 1: (Color online) The total (black), neutrons (blue) and 
protons (red) radial density profiles for asysoft (solid lines) 
and asysuperstifT (dashed lines). 



consider a Skyrme-like [SKM* 



parametrization 
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for q — n{p) and po denotes the 
saturation density. The saturation properties of sym- 
metric nuclear matter are reproduced with the values of 
the coefHcients A = -356MeV, B = 303AfeF, a = 1/6, 
leading to a compressibility modulus K — 200Af eV . For 
the isovector sector we employed three different param- 
eterizations of C{p) with the density: the asysoft, the 
asystiff and asysuperstiff respectively, see [3l| for a de- 
tailed description. The value of the symmetry energy, 
f£ , £iP}_P_ 
3 
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at saturation, as well as the 



slope parameter, L = Spo 



dEsym/A 



are reported in 



dp 

Table |T] for each of these asy-EoS. Just below the satu- 
ration density the asysoft mean field has a weak varia- 
tion with density while the asysuperstiff shows a rapid 
decrease. Then, due to surface contributions to the col- 
lective oscillations, we expect to see some differences in 
the energy position of the dipole response of the system. 

The numerical procedure to integrate the transport 
equations is based on the test-particle (t.p.) method. 
For a good spanning of phase-space we work with 1200 
t.p. per nucleon. We consider the neutron rich nucleus 
^^^Sn and we determine its ground state configuration 
as the equilibrium (static) solution of Eq.([T|). Then pro- 
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ton and neutron densities pq{r,t) 



(27rft)3- 

can be evaluated. The radial density profiles for two 
asy-EOS are reported in Fig. [T] As an additional check 
of our initialization procedure, the neutron and proton 
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mean square radii (r ) 



the skin thickness AR„ 
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Pq{r,t)d r, as well as 



were also cal- 



culated in the ground state and shown in Table HI The 
values obtained with our semi-classical approach are in a 
reasonable agreement with those reported by employing 
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4.65 


0.31 



TABLE I: The symmetry energy at saturation (in MeV), the 
slope parameters, neutron rms radius, protons rms radius, 
neutron skin thickness for the three asy-EoS. 



other models for similar interactions [32!|. The neutron 
skin thickness is increasing with the slope parameter, as 
expected from a faster reduction of the symmetry term 
on the surface [1, |3lj . This feature has been discussed in 
detail in 0. 

To inquire on the collective properties of the pygmy 
dipole we excite the nuclear system at the initial time 
t = to = 30/m/c by boosting along the z direction all ex- 
cess neutrons and in opposite direction all core nucleons, 
while keeping the CM of the nucleus at rest (Pygmy-like 
initial conditions). The excess neutrons were identified as 
the most distant — 32 neutrons from the nucleus CM. 
Then the system is left to evolve and the evolution of the 
collective coordinates Y, Xc and X, associated with the 
different dipole modes, is followed for 600 fm/c by solving 
numerically the equations ([T]). The simple estimate of the 
EWSR provided by the HOSM suggests, when compared 
to experiments, that some of the neutrons boosted in 
the initial conditions are still bound to the core. This 
is confirmed by the transport simulations. Indeed, apart 
from the quite undamped oscillations of the Y coordi- 
nate, we also remark that the core does not remain inert. 
In Fig. [5] we plot the time evolution of the dipole Dy, of 
the total dipole D and core dipole Dc moments, for two 
asy-EoS. As observed, while Dy approaches its maximum 
value, an oscillatory motion of the dipole Dc initiates and 
this response is symmetry energy dependent: larger is the 
slope parameter L, more delayed is the isovector core re- 
action. This can be explained in terms of low-density 
(surface) contributions to the vibration and therefore of 
the density behavior of the symmetry energy below nor- 
mal density: a larger L corresponds to a larger neutron 
presence in the surface and so to a smaller coupling to 
the core protons. We see that the total dipole D{t) is 
strongly affected by the presence of isovector core oscilla- 
tions, mostly related to the isovector part of the effective 
interaction. Indeed, D{t) gets a higher oscillation fre- 
quency with respect to Dy, sensitive to the asy-EOS. The 
fastest vibrations are observed in the asysoft case, which 
gives the largest value of the symmetry energy below sat- 
uration. In correspondence the frequency of the pygmy 
mode seems to be not much affected by the trend of the 
symmetry energy below saturation, see also next Fig. [3l 
clearly showing the different nature, isoscalar-like, of this 
oscillation. For each case we calculate the power spec- 
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FIG. 2: (Color online) The time evolution of the total dipole 
D (top) of the dipole Dy (middle) and of core dipole Dc 
for asysoft (blue, solid) and asysuperstiff (red, dashed) EoS. 
Pygmy-like initial excitation. 
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FIG. 3: (Color online) The power spectrum of total dipole 
(left) and of the dipole Dy (right) (in fm} /(?), for asysoft 
(blue, solid line), asystiff (black, dot-dashed line) and asysu- 
perstiff (red, dashed line) EoS. Pygmy-like initial conditions. 



larly for D. The results are shown in Fig. [31 The position 
of the centroids eorresponding to the GDR shifts toward 
larger values when we move from superasystifF (largest 
slope parameter L) to asysoft EoS. This evidences the 
importance of the volume, S-J component of the GDR in 
^^^Sn. The energy centroid associated with the PDR is 
situated below the GDR peak, at around S.SMeF, quite 
insensitive to the asy-EOS, pointing to an isoscalar-like 
nature of this mode. A similar conclusion was reported 
within a relativistic mean- field approach (ssj . While in 
the schematic HOSM all dipole modes are degenerate, 
with an energy E = 41/1^^/'^ « %MeV for ^^^5*71, within 
the Vlasov approach the GDR energy is pushed up by the 
isovector interaction. Hence the structure of the dipole 
response can be explained in terms of the development 
of isoscalar-like (PDR) and isovector-like (GDR) modes, 
as observed in asymmetric systems [s^l. Both modes 
are excited in the considered, pygmy-like initial condi- 
tions. Looking at the total dipole mode direction, that 
is close to the isovector-like normal mode, one observes 
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FIG. 4: (Color online) The same as in Fig. [3] but for a GDR- 
like initial excitation. 



a quite large contribution in the GDR region. On the 
other hand, considering the Y direction, more closely re- 
lated to the isoscalar-like mode, a larger response am- 
plitude is detected in the pygmy region. To check the 
influence of the initial conditions on the dipole response, 
let us consider the case of a GDR-like excitation, corre- 
sponding to a boost of all neutrons against all protons, 
keeping the CM at rest. The initial collective energy 
corresponds to first GDR excited state, around 15MeV. 
Now the initial excitation favours the isovector-like mode 
and even in the Y direction we observe a sizeable contri- 
bution in the GDR region, see the Fourier spectrum of 
Dy in Fig. 21 From this result it clearly emerges that 
a part of the excess neutrons is involved in a GDR 
type motion and the relative weight depends on the sym- 
metry energy: more neutrons are involved in the pygmy 
mode in the asysuperstiff EOS case, in connection to the 
larger neutron skin size. We have also checked that, if 
the coordinate Y is constructed by taking the Ny most 
distant neutrons (with Ny < Ne), the relative weight in- 
creases in the PDR region. In any case, since part of the 
excess nucleons contributes to the GDR mode, a lower 
EWSR value than the HOSM predictions corresponding 
to Ny = Ne is expected. Indeed, in the Fourier power 
spectrum of D in Fig. 21 a weak response is seen at the 
pygmy frequency. 

These investigations also rise the question of the ap- 
propriate way to excite the PDR. Nuclear rather than 
electromagnetic probes can induce neutron skin excita- 
tions closer to our first class of initial conditions [35| . 
In the case of the GDR-like initial excitation we can re- 
late the strength function to Im{D{uj)) [s^ and then the 
corresponding cross section can be calculated. Our esti- 
mate of the integrated cross section over the PDR region 
represents 2.7% for asysoft, 4.4% for asystiff and 4.5% 
for asysuperstiff, out of the total cross section. Hence 
the EWSR exhausted by the PDR is proportional to the 
skin thickness, in agreement with the results of [s^l- The 
fraction of photon emission probability associated with 
the PDR region can be estimated from the total dipole 
acceleration, within a bremsstrahlung approach [38| . We 



obtain a percentage of 4.7% for asysoft, 7.7% for asystiff 
and 9% for asysuperstiff EOS, consistent with the previ- 
ous interpretation. 
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Summarizing, in this work we evidence, both within 
HOSM and a semi-classical Landau- Vlasov approach, the 
existence, in neutron rich nuclei, of a collective pygmy 
dipole mode determined by the oscillations of some excess 
neutrons against the nuclear core. From the transport 
simulation results the PDR energy centroid for ^^^Sn ap- 
pears around 8.5 MeV, rather insensitive to the density 
dependence of the symmetry energy and well below the 
GDR peak. This supports the isoscalar-like character of 
this collective motion. A complex pattern, involving the 
coupling of the neutron skin with the core dipole mode, 
is noticed. While HOSM can provide some predictions 
of the EWSR fraction exhausted by the pygmy mode, 
N Z 

fy = ^ , depending on the number Ny < Ne of neu- 

trons involved, the transport model indicates that part 
of the excess neutrons are coupled to the GDR mode 
and gives some hints about the number of neutrons, Ny, 
actually partecipating in the pygmy mode. This reduces 



considerably the EWSR acquired by the PDR, our nu- 
merical estimate providing values well below 10%, but 
proportional to the symmetry energy slope parameter L, 
that affects the number of excess neutrons on the nuclear 
surface. We consider these effects as related also to the 
S-J component of the dipole dynamics in medium-heavy 
nuclei. It is therefore interesting to extend the present 
analysis to lighter nuclei, like Ni or Ca isotopes, where 
the Goldhaber- Teller component can be more important. 
We would like to mention that such self-consistent, trans- 
port approaches, can be valuable in exploring the collec- 
tive response of other mesoscopic systems where similar 
normal modes may manifest, see [39|. 
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